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We study the time dependent electron-electron and electron-hole correlations in a mesoscopic 
device which is splitting an incident current of free fermions into two spatially separated particle 
streams. We analyze the appearance of entanglement as manifested in a Bell inequality test and 
discuss its origin in terms of local spin-singlet correlations already present in the initial channel and 
the action of post-selection during the Bell type measurement. The time window over which the 
Bell inequality is violated is determined in the tunneling limit and for the general situation with 
arbitrary transparencies. We compare our results with alternative Bell inequality tests based on 
coincidence probabilities. 



I. INTRODUCTION 

Quantum entanglement of electronic degrees of free- 
dom in mesoscopic devices has attracted a lot of interest 
recently. The first proposal to probe localized entangled 
electrons through transport and noise measurements* 
soon lead to specific structures which generate spatially 
separated streams of entangled particles^. One class 
of devices makes use of a superconducting source emit- 
ting Cooper pairs into a normal-metal structure with two 
leads in a fork geometry: entanglement has its origin 
in the attractive interaction binding the electrons into 
Cooper pairs, while the spatial separation of correlated 
electrons is arranged for by suitable 'filters—. Another 
class of devices makes use of Coulomb interactions in 
confined geometries^. All of these proposals involve elec- 
tronic spins as the entangled quantum degrees of free- 
dom; an alternative scheme has been pointed out by 
Samuelsson et alA who propose a setup where real space 
orbital degrees of freedom become entangled. Besides 
these proposals for the generation of spatially separated 
entangled pairs, the implementation of Bell inequality 
tests probing their entanglement has been discussed in 
detail^*&. The combination of sources for the creation 
and methods for testing the correlations of entangled 
particle streams are first steps towards establishing this 
quantum resource for solid state based quantum infor- 
mation technology. 

An interesting proposal has recently been made by 
Beenakker et ali (see also Ref. and the note|!j): using 
a two-channel quantum Hall device with a beam splitter, 
they suggest a setup generating two streams of entangled 
electron-hole pairs and confirm the presence of correla- 
tions through a Bell inequality (BI) test. A crucial differ- 
ence to previous proposals is the absence of interactions 
generating the entanglement (see also Refs. I^ITolEll) . In 
this paper, we investigate a similar setup involving a 
mesoscopic normal-metal structure arranged in a fork 
geometry and generating two streams of correlated elec- 
trons in the two arms of the fork, see Fig.^ In our formu- 
lation of the Bell inequalities we take special care to use 
only directly measurable observables. We find that a Bell 
inequality test involving correlations between the elec- 



ts v 




t *T o 



source XJlJl 
lead 




(a) 



FIG. 1: (a) Mesoscopic normal-metal structure in a fork ge- 
ometry generating two streams of spin-correlated electrons in 
the two arms of the fork. The Bell type setup detects spin- 
currents 7j, i = 1, 3, projected onto the directions ±a in the 
upper arm and correlates them with spin-currents Ij, j = 2, 4, 
projected onto the directions ±b in the lower arm. (b) In the 
tunneling limit (with a small transparency Td < 1 in the 
'down' arm) a fraction of spin-entangled electrons is split into 
the two arms (components /3 and 7); their correlations can 
be efficiently measured in a Bell type setup sensitive to hole 
(electron) currents in the 'up' (down) lead. 



tronic spin-currents in spatially separated leads exhibits 
violation at short time scales r < r BI , with T m given by 
the single particle correlation time ry = Ti/eV . We show, 
that the time t bi can be considerably extended in the tun- 
neling limit, where the propagation into one of the arms 
(we choose the 'down' lead in Fig.Q} is strongly reduced 
by a small transparency This small transparency can 
be exploited by a suitable choice of the observables in the 
Bell inequality test: going over to hole currents in the 
well conducting 'up' lead, we find the long violation time 
t B i = T~v/\/Td- We also analyze an alternative Bell in- 
equality based on coincidence probabilities derived from 
electron-hole number correlators and find it violated on 
even longer times r BI = tv/T^. When tracing the origin 
of these BI violations, we find that the fermionic statis- 
tics already enforces the injection of spin-singlet corre- 
lated electron pairs with correlations extending over the 
distance i> F ry with v F the Fermi velocity. In addition, the 
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splitting of this pair and a subsequent post-selectioniiiiS 
in the Bell type measurement are crucial steps in ren- 
dering the original spin-entanglement amenable to obser- 
vation and potentially useful as a resource of entangled 
quantum degrees of freedom. 

In the following, we first analyze (Sec. II) the two- 
particle density matrix of the injected particle stream 
and find it to be locally singlet correlated; the extension 
to the scattered state behind the beam splitter demon- 
strates that this entanglement is preserved. We proceed 
with a discussion of Bell inequality measurements (Sec. 
Ill) and compare two types of Bell parameters, one based 
on current cross-correlators and a second one starting 
from coincidence probabilities, i.e., particle number cor- 
relators. We then turn to the actual calculation of the 
finite time spin-current correlators and their combination 
into the Bell parameter (Sec. IV). Results are presented 
for Bell inequalities involving electron-electron (general 
case) and electron-hole (tunneling limit) current correla- 
tions as well as for the corresponding expressions based 
on coincidence probabilities; we demonstrate that the vi- 
olation of the Bell inequality depends sensitively on the 
choice of observables. Finally, we present our conclusions 
in Sec. V. 



II. TWO-PARTICLE CORRELATIONS 

In order to analyze the properties of the injected elec- 
trons, we determine the two-particle density matrix (or 
pair correlation function) within the source lead, 

g ff (x, y) = Tr(p ¥ ai (x)¥ a2 (y)i> a3 (y)K 4 (*)) (1) 
= Wi\a 4 )(a 2 \a 3 ) G 2 (0) - (<7i|<7 3 )<ff 2 |ff4>|G(a; - y)\\ 

where we have introduced the one-particle orbital corre- 
lator G(x) = (# t (a;)>I'(0)), with $(x) the field operator 
describing spinless electrons in the source lead and o~i are 
spin indices. We ignore contributions from backscatter- 
ing originating from the splitter and describe the source 
lead as a ballistic wire connecting two reservoirs with 
Fermi levels shifted by the voltage bias eV. The correla- 
tor can be separated into an equilibrium- and an excess 
part vanishing at zero voltage, G(x) = G eq (x) + G ex (x), 

G oq (x) = [sin(fc F x)]/7ra, (2) 
G cx {x) = e -*(fe+fev/2)x [ s i n { kv x/2)]/TTX, (3) 

with k v — k F (eV/2e F ). The equilibrium part of the 
pair correlator g c ff q (x,y) = (CTi|cr 4 )(cr 2 |cr3)[G oq (0)] 2 - 
(ci|cr3)((72|c4)|G eq (j; — y)\ 2 describes the exchange cor- 
relations between two fermions with singlet correlations 
decaying on the Fermi wave length A F . The excess part 

g?(x,y) = H(7 4 ) (a 2 |a 3 )[G cx (0)] 2 

-{a l \a,){a 2 \a i )\G^{x-y)\ 2 (4) 

associated with the additional injected electrons ex- 
hibits singlet correlations on the much larger scale X v = 



A F (e F /eV) (an additional mixed term describing the de- 
formation of the equilibrium exchange holei^ due to 
the bias is not relevant to our discussion and we ig- 
nore it here). We conclude, that the excess particles in- 
jected by the reservoir form a stream of singlet- correlated 
pairs. Furthermore, again due to the Fermi statistics, 
these singlet-pairs propagate as a regular sequence of 
wave packets separated by the single particle correla- 
tion timei^ Ty = h/eV corresponding to the singlet cor- 
relation length A v . The singlet-pairs are conveniently 
described through the state \^f; n ) = l e |j)s (with in- 
dex 's' for 'source') or the corresponding wave func- 
tion Vfin = 4>l(j)lxll: with the spin-singlet wave function 
Xsg = \x)x\ ~ XjXtIA/2 and upper indices 1 and 2 de- 
noting the particle number. Note that the emission of 
singlet-pairs from the normal reservoir follows naturally 
from the identical orbital wave function describing elec- 
trons with opposite spins within the reservoir. 

Next, we analyse the scattered state propagating in the 
two leads of the fork in Fig. ^ The pair correlation func- 
tion describing particles propagating in different leads (x 
in 'u' and y in 'd') takes the form 

g 3 (x,y) = (ai|a 4 )(<72k 3 )G u (0)G d (0) 

-{<yiW 3 ){a 2 \a i )\G ud {x-y)\\ (5) 

and has to be calculated with the scattering states 




where a e , 6 e , and c e denote the annihilation operators 
for spinless electrons at energy e in leads 'u', 'd', and 
's' and with a time evolution oc exp(— iet/h) (here, t su 
(fdu) an d i s d (tud) describe particle transmission from 
the source (down) lead into the 'up' lead and from the 
source (up) lead into the 'down' lead; r u , r d denote the 
reflection amplitudes into leads 'u' and 'd'). The or- 
bital one-particle correlators G x (x) = (^(^^(O)) = 
G eq (x) + TxG ex (x) describe particles with coordinate x 
residing in the same lead 'x' equal to 'u' or 'd' (T u = \t su \ 2 
and Td = |t s d| 2 describe the transmission probabili- 
ties into the 'up' and 'down' lead). The one-particle 
cross correlator between leads 'u' and 'd' takes the form 
G ud (a; -y) = (*£(a;)*d(lO> = t* su t sd G e *(x - y) with co- 
ordinates x and y residing in the leads 'u' and 'd'. The 
excess part of the pair correlation function reads 

g?(x,y) = r u ^ d [<ctx | CT4 ) <ct 2 J C7 3 > [G ex (0)] 2 

-<aik3><^k4>|G e *(z-y)| 2 ]. (6) 

This result is identical in form with the excess pair corre- 
lator 10} in the source lead, however, it now describes the 
correlation between a singlet-pair split into the leads 'u' 
and 'd'. Note the preservation of the singlet-correlations 
which are maximal for x = y and decay on a distance 
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| a; — y\ ~ X v , where the coordinates x and y belong 
to different leads. The scattering state describing the 
propagation of the singlet-pair behind the splitter can be 
written in the form 

mit) =4k4l)ul°)d+^d|0)u|e s T s i ) d (7) 

+V / 2t S u4d[|et)u|e;)d - |ej,) u |e T ) d ] , 

where the first two terms describe the propagation of the 
singlet-pair |e^) x in leads 'x' equal 'u' or 'd', with a wave 
function of the form ^^xXsg- The last term describes a 
singlet-pair split between the 'up' and 'down' leads with a 
wave function [0„<^> d + 4'd ( t ) u}xll ■ A coincidence measure- 
ment of electrons in leads 'u' and 'd' projects the scat- 
tered state l^out) onto this spin-entangled component 
with spatially separated electrons in leads 'u' and 'd'. 

In the tunneling limit 7 (T u ~ 1 and T d <C 1) most of 
the incoming singlet pairs propagate into the well con- 
ducting 'up' lead and only rarely (with amplitude i su i s d) 
split into both leads. The absence of an electron in the 
'up' lead then manifests itself as the presence of a hole 
and it is favorable to go over to a hole representation, 

,hc \ _ 



Ie s ^) d 



,,:l"' T x/u|e Ti /d (8) 



t 2 jo) u \o) d + ejh 

v / 2i su i s d[|^i)u|ex)d - | /it) u| e T/d 



and the hole current = (2e 2 /h)V — j„. The first term 
in (JSJ) (the component a in Fig. 1(b)) describes a filled 
Fermi see in the upper lead combined with a vacuum 
state in the lower lead, and hence no particle can be 
detected. The second term (component 6) accounts for 
the rare processes where both electrons propagate to the 
'down' lead; its contribution spoils the maximum vio- 
lation of the BI and restricts the use of the tunneling 
limit. The most relevant terms are the last two (0 and 
7) describing the splitting of the singlet-pair between 
the two leads and the formation of a spatially separated 
spin-entangled electron-hole pair with the hole (electron) 
propagating in the upper (lower) lead; this electron-hole 
component is detectable in a coincidence measurement 
using a hole (particle) detector in the upper (lower) lead. 



III. BELL INEQUALITIES 

A. Bell inequality with current correlators 

The original goal in setting up the Bell inequalities^ 
was to devise a scheme allowing for the differentiation be- 
tween classical correlations appearing in a local hidden 
variable theory and non-local correlations as they show 
up within a quantum mechanical framework. Accord- 
ingly, early experiments in optics addressed these funda- 
mental questions dealing with the validity of quantum 
mechanics. Recently, Bell inequalities have been dis- 
cussed in the context of mesoscopic electronic devices. 
One should admit, that the corresponding electronic se- 
tups are probably less suitable for addressing fundamen- 
tal issues of quantum mechanics. In a more pragmatic 



approach, the Bell inequalities in mesoscopic physics are 
used as a test for the presence of entanglement or even 
for a quantitative measurement of the degree of entan- 
glement between quantum degrees of freedom. 

Defining appropriate Bell inequalities in mesoscopic 
physics is nevertheless a non-trivial issue as those ob- 
servables suitable for direct measurement in optics are 
not necessarily available in mesoscopics; this is why we 
attempt an extended and detailed analysis below. In this 
context, we keep the discussion on a level where funda- 
mental and practical issues can be easily identified. 

The explicit form of the Bell inequality we are going to 
use in the present paper has been introduced by Clauser 
and Horna^ based on the original discussion of Bell 15 . 
It derives from the lemma saying that, given a set of real 
numbers x, x, y, y, X, Y with |x/X|, and 

\y/y\ restricted to the interval [0, 1], the inequality 



\xy — xy + xy + xy\ < 2\XY\ 



(9) 



holds true. In the Bell type setup of Fig. ^a) one 
measures the spin-projected electronic currents If, i = 



1,. 
If 



, 4, and defines the quantities x = Jf 



it x 



Jf and y = Jf 



12, Y 



I? for fixed orienta- 



tions a and b of the polarizers (and similarly x and y for 
orientations a and b). Our Bell setup then measures the 
correlations 



Cr;(a,b;r) = {If(r)I](0)) x 



(10) 



J dAp(A)/f(a,A;r)/ J e (b,A;0) 



between the spin-currents If 



1,3, in lead 'u' pro- 



jected onto the directions ±a and their partners /?, 
j = 2,4, in lead 'd' projected onto ±b. Within a local 
hidden variable (A) theory, the average (. . .)\ is taken 
with respect to the distribution function p(A); specify- 
ing a theoretical framework such as quantum mechanics, 
these averages are replaced by quantum mechanical av- 
erages. Using the above definitions of x, y, X, and Y, we 
obtain the current difference correlator 



£(a,b;r) - 



([i?(T)+IS(r)][JS(0)+JS(0)]) A : 
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and evaluating these for different combinations of direc- 
tions a, a, b, and b we can combine them into the Bell 
inequality 

£ BI (T) = |£(a,b)-£(a,b)+£(a,b)+£(a,b)| < 2. (12) 

We can further process the difference correlator Hll|) 
and separate the current correlators into an irreducible 
part C°°(a,b;r) = (<5I?(r) 8If(0)) x with SIf{r) = 
If(r) — (If)\ and a product of average currents and 
rewrite E(a, b; r) in the form 



£(a,b;r) 



Ccc pec pec 1 pee i A ec 
12 °14 — °32 ' °34 T yv - 



Pee 1 pec _i_ p< 

°12 "+" °14 ' °; 



32 



G 



34 



A^_ e ' 



(13) 
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withA c ± c = [</?) A ±</f) A ][(/f) A ±(/|) A ]. The average cur- 
rents are related via (If}\ = (-Tf) A = (^u)a/2 and (/|) A = 
(^4>a = (IZ)x/2 and thus Ai c = 0, A- = (I°)x(I c d )x, 



Ccc /^cc /"fee j_ /"^ec 
12 °14 ~ °32 ~i~ °34 

Cfi + Cfi + Cfi + Cfi + (/s)a(J3>a • 



(14) 



In the tunneling limit, the electronic currents i| in the 
'up' lead are replaced by hole currents if — I^ ax — 
If, where i^ax denotes the current in the open chan- 
nel (within a quantum mechanical framework, I^ ax = 
(e 2 /h)V). 



B. Bell inequality with coincidence probabilities 

The Bell inequality (|12Jl explicitely depends on the 
delay time r which appears naturally through the time 
dependence in the current correlators in (|llf> and 1)14)1. 
This implies that the violation of the Bell inequality de- 
pends on the delay- or measurement time of the exper- 
iment, a feature not encountered in traditional optical 
setups. In optics, the quantity usually involved in this 
type of analysis is the coincidence probability (a, b) 
for the simultaneous detection of two photons with po- 
larizations along a and b. The indices i and j specify 
the directions along and perpendicular to the polarizer 
— here, these correspond to the spin-up and spin-down 
events triggering a signal in the detectors i — 1,3 or 
j = 2,4. In optics, these coincidence probabilities are 
directly measurable and can be combined into a Bell in- 
equality E BI = \E(a,b)-E(a,b)+E(a,b)+E(a,b)\ < 2 
with 

E{a,b) = F 12 (a,b) - P u (a,b) - P 32 (a, b) + P 34 (a, b). 

Contrary to the situation in optics (where photons are 
annihilated in the detector), the coincidence probability 
Pij (a, b) is not a directly measurable quantity in meso- 
scopics, where the observables of choice are the (elec- 
tron or hole) currents I z e ' h . The expression of Pij (a, b) 
through measurable currents then requires some care and 
we provide a detailed discussion here in order to compare 
the approach based on coincidence probabilities with the 
one based on current correlators, see section III. A. 

A natural way to define a coincidence probability in 
mesoscopics is through particle number correlators 



^(a,b;r) = <JV?(r)J^(T))A 



where 



N!(r) 



dtim 



is the number of electrons counted by the detector i dur- 
ing the accumulation time r. Here, we are interested in 
the simultaneous detection of two particles, one appear- 
ing in the upper lead 'u', the other in the down lead 'd'. 
In order to obtain such a coincidence probability from the 



number correlator i|15[) we have to restrict the accumula- 
tion time such that only events {0, 0} (no particles) {1,0} 
(one particle in detector i), {0, 1} (one particle in detec- 
tor j), and {1, 1} (two particles, one in detector i and one 
in j) are accounted for; out of these, the coincident events 
{1, 1} then contribute to ACJJ. Events of the type {2, 1}, 
{1,2}, ... have to be avoided through proper time limi- 
tation. Restricting the accumulation time r to a value T\ 
such that at most one particle is counted, (N? (ti)) a < 1, 
and using a proper normalization, we can find an ex- 
pression for the coincidence probability in terms of the 
number correlator 1151) . 



i^(a,b;ri) 



ac? 



AC? 



ACS 



AC 



(17) 
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with all correlators AC|J(a, b, t\) evaluated at a fixed time 
Ti and fixed directions a and b, cf. (|15|) . The condition 
(N£(ti)) < 1 requires that T\ is smaller than the time 
tj = e / maxi[(If) \] between subsequent events, n <C Tj. 
Note, that during the short accumulation time n less 
than one particle contributes on average, however, this 
is corrected for by the proper normalization. Also, we 
note that our definition (|17fl for the coincidence rate is 
not i dentic al to the one introduced in Ref. and used in 
Refs. l7illt we will return to this point later. 

The quantity l|17f) can be further processed in the limit 
of rare events, i.e., in the tunneling limit (Tj <C 1). In 
this situation it is advantageous to go over to electron- 
hole correlators^ with the hole current defined via if — 
^max — It (note that in Ref. |6| the electronic currents in 
the two normal leads originate from the low-rate emission 
of Cooper pairs from the superconductor and hence are 
both small; alternatively, one may view these excess cur- 
rents as arising from Andreev-scattering at the normal- 
superconductor interface where an incident electron is 
reflected as a hole). Again, the electron-hole number 
correlators entering (|17[) are split into an irreducible and 
a reducible part, 

K^{T 1 )=K^{r)+{lf) x {It) x rl. (18) 

with the irreducible number correlator defined as 

K${r) = (5NHt)5N?(t))x = fait fdt' C%(t>) (19) 

Jo J-t 

Next, we can make use of the fact that the irreducible 

(15) current correlator Cf^(i') rapidly decays on a time scale 
t c (within a quantum mechanical framework, t c = ry = 
fi/eV is the single particle correlation time). For t\ ^S> 
t c , we can approximate the irreducible part of AC^ e (ri) 

(16) by the time-independent zero-frequency noise correlator 



Sf?(u = 0) (cf. Ref. 



4f(n) = Sfnu; 



0)n + (7^) A (7| 



(20) 



Inserting this expression into the correlators E we ar- 
rive at the expression analogous to the one introduced in 
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In a final step, we demonstrate that we can ignore 
the current product term A^n in the tunneling limit. 
In order to gain more insight, we analyze the situa- 
tion theoretically within a quantum mechanical frame- 
work for non-interacting electrons: with Sf°(ui = 0) = 
\( ai \b 3 )\ 2 (e 2 /h)T u T d eV and (if) a (Ip = (e 2 /h)T d V 
(we assume that 1 — T u ~ T d ; the vectors and b, de- 
note the directions associated with the detectors i and 
j) we find that K%( n ) = {e 2 /h)nT d eV(\(a i \b j )\ 2 T u + 
TiXd/27rry). The second term is due to the current prod- 
uct term and can be dropped provided that T\T d j2-KTy <C 
1 (the appearance of the small factor T d is due to the use 
of hole currents). In the tunneling limit, Ty/T d = 77 3> 
Ty and we find a large time window t c — Ty <C t\ <C 77 
within which we can choose our accumulation time t\ 
such that both the definition l|17(l and the separation l|20() 
can be properly installed and the current product term 
is small. It is crucial to understand that under these cir- 
cumstances the accumulation time 77 appears as a the- 
oretical quantity which is needed only in the derivation 
of the Bell parameter; once we have demonstrated that 
77 can be chosen such that the current product term is 
small, the latter can be dropped and 77 disappears from 
the Bell parameter. 

While the above idealized theoretical consideration 
serves as a guideline, the situation in a real experiment 
may be complicated due to interactions and other effects 
which change the value of the correlation time r c , i.e., r c 
may differ substantially from Ty. Therefore, in an actual 
experiment testing the Bell inequalities the decay time 
r c of the correlator should be measured independently in 
order to verify that the second term in (|20l) involving the 
product of average currents is indeed small on this time 
scale and thus can be dropped. Fortunately, in the tun- 
neling limit the admissible time window r c <C 77 <C 77 
is large and this test can be carried out with only an 
approximate knowledge of the correlation time t c . 

Once we have verified that we can drop the correc- 
tion term (I^)\(Ii)\T 2 in (|2*U)l we can express the co- 
incidence probabilities through the (time-independent) 
zero- frequency noise correlators Sf°(ui — 0) alone, 
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the correlator E(a, b) then takes the simple form 
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The results J22l> and now are independent of time, 
although the original expression l|17[) involved the time 
restriction n <C 77; their use is restricted to the tunneling 



limit as only in this case the correction term ex A^ e in the 
denominator of E can be dropped. 

Let us next concentrate on the general situation away 
from the tunneling limit which turns out quite different. 
In this case, there is no advantage in going over to hole 
currents and we work with electronic currents in both 
leads. Second, in the absence of a small tunneling prob- 
ability we have no separation of time scales and r c ~ 77 
are of the same order, hence the accumulation time T\ 
has to be smaller than (or at most of the order of) the 
single particle correlation time r c , t\ < r c . As a conse- 
quence, one cannot express the coincidence probability 
P^ e through the zero- frequency noise correlator Sfj. In- 
stead, we make use of the irreducible number correlator 
Kfj(r), cf. (|19f) . and rewrite the coincidence probability 
in the form 



fS e (a,b;n) = 



Kf 
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(24) 



Combining the coincidence probabilities into the expres- 
sion E(sl, b), the products of average currents cancel in 
the numerator, Aff = 0, however, these products do not 
cancel in the denominator and restrict the violation of 
the Bell inequality to small time intervals (see below). 
The final correlator entering the Bell parameter E B1 then 
takes the form 



£(a,b;ri) 



K ee — K cc 
iv 12 j\ 14 



K 



32 



IC 



34 



K\l + Kfl + KH + Kit + (I°)x(IZ)xt 2 



(25) 

Comparing l|14fl with l|25|l it turns out that both Bell pa- 
rameters E B1 and E BI , when derived carefully, are based 
on the same correlations, once expressed directly through 
currents, the other time through number correlators. In 
particular, both Bell parameters depend on time and, as 
we will see below, the violation of the Bell inequalities is 
limited to times n < t c ~ 77. Quite interestingly, this 
coincides with the time restriction imposed by our defini- 
tion l|17l) of the coincidence probability in terms of num- 
ber correlators — it turns out that the Bell inequality can 
be violated during those times T\ where the (normalized) 
number correlator assumes the meaning of a coincidence 
probability. 

The correlator 125fl has first been introduced and used 
in Ref.^in order to analyze the entanglement of electrons 
injected by a superconductor into a normal metal fork. 
On the other hand, the simplified version (|23|l valid in 
the tunneling limit and first introduced (via an alterna- 
tive route) in Ref. has been successfully applied in Ref. 
7. However, its application in Ref. to a situation away 
from the tunneling regime led to the appearance of a 'spu- 
rious amplification factor' in the relation between the Bell 
parameter and the concurrence, a consequence of ignor- 
ing the presence of the current product term oc (/u)(^d)- 
Instead, the subsequent use& of the full expression (f5Sf> 
removed this problem and established the violation of 
electronic Bell inequalities in the non-tunneling regime 
at short times r < Ty — Ti/eV. 
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In the discussion above, we have been careful to sepa- 
rate theoretical from experimental input in the construc- 
tion of the Bell parameter and have provided final ex- 
pressions involving only experimental input. Depending 
on the chosen variables and on the physical situation, we 
have seen that a short time measurement is required in 
general; in the tunneling limit, an approximate determi- 
nation of the correlation time r c together with the zero- 
frequency noise correlator is sufficient. It is interesting to 
compare our point of view expressed in the above deriva- 
tion with the approach introduced in Ref. 0, where the 
authors derive a Bell parameter which does not involve 
a short time measurement. While this scheme works the 
same way as ours in the tunneling limit (rough estimate 
of t c and knowledge of the zero-frequency noise are suf- 
ficient), it does require the additional precise knowledge 
of the correlation time r c away from the tunneling limit, 
which cither requires an accurate theoretical evaluation 
for the device at hand or again necessitates a short time 
measurement. 

In Refs. IHlH the coincidence probability Pjj(a, b) has 
been defined as an equal time correlator of the form, 

Py(a,b) cx (*)&(*)) (26) 

ex = (if)(i°) + (si?(t)si°(t)). 

Here, the operators a* and bj annihilate electrons in the 
upper and lower leads in the detectors i = 1,3 and 
j = 2, 4. In contrast to our definition 117|) above, which 
has been based on physically measurable quantities and 
which involves a time restriction T\ -C tj, the definition 
((26(1 is a time independent correlator but probably can- 
not be measured in a mesoscopic setting. The idea put 
forward in Ref. ^3 then is, to use quantum mechanics to 
reexpress this theoretical definition through measurable 
quantities. In the tunneling limit, after transformation to 
electron-hole currents, the current product term in 1(26|) 
can be dropped (after proper experimental check, see the 
discussion above) and the irreducible current correlator 
can be expressed through the zero-frequency noise corre- 
lator via 



made in Ref. I 111 then is to construct the coincidence prob- 
ability from the combination 



with the correlation time r r defined via 



= / dt 



'(«h(0)<5/ c (0))' 



(27) 



(28) 



Assuming that r c does not depend on the lead indices 
i and j, this time scale disappears after normalization 
and one arrives at the formula 1)23(1 expressed through 
the zero frequency noise Sff as derived above starting 
from the measurable expression l(17|l for the coincidence 
probability. 

However, away from the tunneling limit, the current 
product term in 1(26(1 cannot be dropped and the corre- 
lation time t c does not vanish any longer; the proposal 



Pij oc t c x Sff 



WW 



(29) 



which requires the measurement of the average currents 
(If) and (I?) in addition to the zero-frequency noise cor- 
relator. The problematic step in this construction is the 
need for the precise quantitative knowledge of the corre- 
lation time t c , since this parameter now is part of the 
evaluation of the coincidence probability itself (this is 
different from the above discussion of the tunneling limit 
where a rough knowledge of r c was sufficient in order 
to verify that the current product term can be dropped 
from 1(20(1 ). The appearance of this additional parameter 
is a consequence of expressing the equal time correla- 
tor (Sif(t)SI](t)) = j{duj/27r)Sff(uj), involving all fre- 
quencies, through the zero-frequency value Sfj ^w = 0) 



alone. The point of view put forward in Ref. [llj then 
is that the parameter r c shall be obtained from a theo- 
retical calculation, e.g., r c = Anry in a non- interacting 
system where the irreducible current correlator decays 
oc sin 2 (eVt/2%) / (eVt/2h) 2 (note that in Ref. Ej] this re- 
sult is used in the expression for the coincidence proba- 
bility Pij). The proposal to replace r c in ((29(1 through 
a theoretically calculated quantity avoids the need for a 
short time measurement; on the other hand, one has to 
accept that the coincidence probability obtained in this 
manner is subject to the approximations (such as neglect- 
ing effects of interactions, resonances, etc. present in the 
real experiment) made in the theoretical evaluation of 
t c . Alternatively, one might want to obtain r c , cf. 12811 . 
directly from an experiment; however, as r c now is used 
in the evaluation of the coincidence probability, a pre- 
cise knowledge of this parameter is required and hence 
an accurate measurement of the current correlator has 
to be performed. This then boils down to a short time 
measurement and nothing can be gained. 

Below, we take the point of view that the Bell inequal- 
ities should be built from physically measurable quan- 
tities. Starting from the expression l(14|). we proceed 
with its theoretical evaluation in order to predict the ex- 
pected outcome of such a Bell inequality test within a 
quantum mechanical frame. We first determine the fi- 
nite time current cross-correlator between leads 'u' and 
'd' for a stream of spinless fermions; the generalization to 
the spin-1/2 case is straightforward. We express the BI 
in terms of these finite time correlators and find its vi- 
olation for the general case expressed through electronic 
correlators and for the tunneling case involving electron- 
hole correlators. We also derive the results expected from 
the alternative formulation based on coincidence proba- 
bilities. 
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IV. CURRENT CROSS-CORRELATORS 
A. Bell inequalities with electron currents 

Starting from the field operators and de- 
scribing the scattering states in the leads, we deter- 
mine the (electronic) irreducible current cross correlator 
C% b v {t) = (SI°(T,x)SI°(0,y)} with positions x and y in 
the leads 'u' and 'd' using the standard scattering theory 
of noiseii and split the result into an equilibrium com- 
ponent C^ y (r; V = 0) and an excess part C| x y (r; V) 



@x,v ( T ) 



e 2 T du 
h 2 

4e 2 T u T d . 2 eV(r-T_) 



h 2 



a{T + T + ,8)+a{T-T+,8)}, (30) 
a(r - T-,6), 



21) 



with a(r,9) = tt 2 9 2 / sinh 2 [nOr/h], t± = (x±y)/v F , and 
9 the temperature of the electronic reservoirs. In order to 
arrive at the result l|5U|l we have dropped termsi£ small 
in the parameter |e' — e|/e F and have used the standard 
reparametrization of the scattering matrix for a three- 
terminal splitter (see Lesovik et al. in Ref. 0). 

Extending the above results to spin-1/2 particles, 
we introduce the spin-projected field operators ^1,3 = 

^3 (T (±a|o')^'ucr (and similar for the 'd' lead). The corre- 
lators C°°(a, b) relate to the result (|3Tj|) for spinless par- 
ticles via C?/(a,b;r) = ^h^Cf^r) with a 1)3 = ±a 
and b2,4 = ±b. The spin-projections derive from the an- 
gle 6* a b between the directions a and b via (±a| ± b) = 
cos 2 (6»ab/2) and (±a| =pb) = sin 2 (6> ab /2) and the BI JEJ) 
assumes the form 



C l% ( r ) [ C0S #ab - COS 6» a6 + COS 6> ab + COS 6> ab ] 



<1. (31) 



Its maximum violation is obtained for the standard orien- 
tations of the detector polarizations, # a b = #ab = # a b = 
7r /4) ^ab = 37r/4 and the BI reduces to 



2C- (r) + </S)<^> 



< 



1 



V2' 



(32) 



In the limit of low temperatures 9 < eV and for large 
distances x = y Tyv ¥ (allowing us to neglect the equi- 
librium part in the correlator C^ y ) the above expression 
(|3*2"|) reduces to the particularly simple form 



2sin 2 (eVY/2?i) 



< 



-75. ( 33 ) 



T 2 (eV/h) 2 -2sm 2 (eVT/2h) ~ ^ 

where we have used that (7°) = (2e 2 /h)T n V and (1%) = 
(2e 2 /h)T d V. We observe that in this limit the Bell in- 
equality is i) violated at short times r < t bi = ry, see 
Fig. 2, ii) this violation is independent of the transparen- 
cies T u , T d and hence universal, and Hi) the product of 
average currents A+ is the largest term in the denomi- 
nator of H31fl and hence always relevant. Note that the 



important quantity appearing in (|32|l is the space and 
time-dependent correlator C° c y (r). The small quantity 
required for the violation of the BI then is the shifted 
time t — t~ < Ty; placing the detectors a finite distance 
apart one may make use of the additional time delay, 
although the time resolution remains unchanged. 

In the low frequency analysis of Ref. |y no violation in 
the BI had been found for a normal injector, in agree- 
ment with the results found here. On the other hand, 
it has been realized in Ref. that a short time mea- 
surement on the scale Ty can exhibit entanglement in a 
normal system. In particular, the proper use of the Bell 
parameter _E BI with E given by (|25|l provided such an 
entanglement away from the tunneling limit, while the 
use of the expression l)23|) led to a 'spurious amplification 
factor' in the relation between the Bell parameter and 
the concurrence. 




eVxIh 



FIG. 2: Bell inequality test for electron-electron (left) and 
electron-hole (right) currents. The thin line at l/y/2 marks 
the critical value above which the Bell inequality is violated. 



Bell inequalities with electron-hole currents; 
tunneling limit 



Next, we consider the tunneling limit and determine 
the outcome of a Bell measurement involving a hole cur- 
rent = (2e 2 /h) V— 1° in the 'up' lead and the electronic 
current j| in the 'down' lead. The cross-measurement in 
different leads implies that the setup is sensitive only to 
the split-pair part of the scattering wave function I^q^) 
which is fully spin-entangled and hence the Bell inequal- 
ity can be maximally violated. In the tunneling limit, 
the correction term cx t 2 d also contributes a signal and 
spoils the maximal violation, ultimately limiting the use 
of the tunneling limit to devices with large enough trans- 
parency T u in the well conducting lead (small enough 
transparency T d in the blocked lead). 

The calculation proceeds as above but now involves the 
electron-hole correlator At) = (51^ (t, x)<5/ d (0, y)) = 
—C° c v (t) and the product of the electron and hole cur- 
rents A\ c = = (2e 2 /h) 2 T d (l - T U )V 2 (again, 
A 1 ! _ rp^g me q ua iity corresponding to l|33|l now 
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reads 



he 



2sin 2 (eVY/2?i) 



< 



2sm 2 (eVT/2h)+Tn 1 {l ~ T u )T 2 (eV/h) 2 ~ V2 



and an illustration of this result is given in Fig. [21 The 
violation of the Bell inequality in the tunneling limit 
exhibits a much richer structure: i) the violation re- 
quires a minimum transparency T u in the upper lead: 
evaluating E^f at r = + , we obtain the condition 
T u > T min = 2/{V2 + 1) w 0.83. u) For T u > T min 
the Bell inequality is violated during times r < r BI = 
v 7 2 - \/2 V T u ( 1 - T u ) - 1 w tv/a/Th, where we have 



Tu)- 1 « n//V7^, 
assumed 1 — T u w Td as is the case for a splitter with 
a small back reflection r s <C 1; this result is different 
from the time limitation noted in Ref. 0. Hi) The BI 
remains un-violated in narrow intermediate regions sep- 
arated by the single particle correlation time ry and de- 
creases slowly oc t -2 with increasing time, iv) The prod- 
uct of average currents Aiff gives a small correction to 
the denominator in Ef^ at short times. The comparison 
with the electronic result is quite striking: the time inter- 
val over which the Bell inequality is violated is extended 
by a factor l/-y/7d ^> 1 and the universality (i.e., the 
independence on the transmissions T u and Td) is lost. 



C. Bell inequalities with number correlators 

Finally, we quote the results obtained for the Bell mea- 
surement based on coincidence pr obab ilities or number 
correlators, cf. Eq. (|25J) and Refs. IfilSt Again, we split 
the number correlator Kf° (cf. Q19[l. we use the elec- 
tronic version) into an orbital- and a spin component, 
Kfj = |(ai|bj)| 2 K° e y {r) and assume the standard set of 
directions a, b, a, and b in order to arrive at the inequal- 
ity 



2ff"(r)+T*<I»><I5) 



< 



1 

71" 



(34) 



Again, this electronic Bell inequality is universal (cf. E< 
(|30|I % 1 and violated at short times r < jy, cf. also Ref. 

The tunneling limit involving the electron-hole num- 
ber correlator K^ e y (r) = (SN^(x,t)5N^(u,t)) is more 
interesting: The Bell inequality takes the form (|34|l but 
with K cc and (J u )(/|) replaced by K hc and (T^K-^)- 
Its evaluation at r = + provides the same condition 
?u > Tmin as found previously for the violation of the 
Bell inequality. In the tunneling limit (1-T U )/T U < 1, 
the number correlator can be estimated at large times 
t > t v and we find K% y (r) « (e 2 /2ir)T u T d (eVT/h); 
the Bell inequality then is violated for even larger times 
t < t bi = t v (V2- 1)ttT„(1 - Tu)" 1 w Ty/T d . Hence we 
see that the violation appears just over those time scales 
t < tj where the number correlator /C^?(t) assumes the 
meaning of a coincidence probability. Note that the time 
dependence found here is lost once we drop the current 



product term r 2 (T u 1 )(/^), taking us to the time indepen- 
dent result corresponding to l|23|) . 



D. Origin of entanglement 

An interesting question concerns the origin of the en- 
tanglement detected in the Bell inequality measurement 
described in the present paper. We note that in Refs. 

the entanglement had been attributed to the elastic 
scattering in the Fermi sea, although the proper selec- 
tion of a projected wave function component during the 
calculation corresponds to a post-selection. In Ref. ^3 
post-selection was noted to be the origin of entanglement; 
such post-selection creating entanglement is a well known 
mechanism in optics 12 . In both of the above cases, the en- 
tangled degrees of freedom originated from independent 
reservoirs. The situation is slightly more complicated in 
the present case: As shown above, our setup involves a 
simple normal reservoir injecting spin-singlet correlated 
pairs of electrons into the source lead which are conve- 
niently described by the wave function ^>l^ = (pl^lxlg- 
These local spin-singlet pairs are subsequently separated 
in space by a beam splitter and detected in a coinci- 
dence measurement. The measurement is only sensitive 
to pairs of particles propagating in different arms, imply- 
ing a post-selection or projection of the scattering wave 
function during which only its cross term describing a 
split spin-singlet pair survives. In this context it is inter- 
esting to note that the incoming local spin-singlet, being 
a simple Slater determinant, is not entangled according 
to the definition given by Schlicmann et ali& However, 
after the beam splitter the orbital wave function </> s is de- 
localized between the two leads, <p s — ► $ = i S u</>u +i s d0d- 
While the scattered state remains a Slater determinant 
^out = < J )1< I )2 Xsgi the singlet correlations now can be ob- 
served in a coincidence measurement testing the cross- 
correlations between the leads 'u' and 'd'. Hence the 
original spin-entanglement is produced by the reservoir, 
but its observation requires proper projection. It is then 
difficult to trace a unique origin for the entanglement 
manifested in the present Bell inequality test. An appro- 
priate setup addressing this question should involve in- 
dependent reservoirs injecting the particles carrying the 
degrees of freedom to be entangled, e.g., particles with 
opposite spin residing in a Slater determinant of the form 
^tn = [^sT^s| — ^si^stl/v^ which is not entangled in the 
spin variable. Such an analysis has been presented in 
Ref. Il9l with the result, that the orbital projection in the 
coincidence measurement is sufficient to produce a spin- 
entangled state. 



V. CONCLUSION 

In conclusion, we find that spin-entangled pairs of elec- 
trons exist and leave their trace in the violation of Bell 
inequalities in a mesoscopic setup even in the absence 
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of interaction. The source of entanglement is traced 
back to the nature of injected electrons forming a regu- 
lar stream of locally singlet-correlated particles combined 
with a post-selectionii*i£ during the Bell type measure- 
ment. The splitter itself does not contribute to the en- 
tanglement of the pair, but fulfills the crucial task of sep- 
arating the spin-entangled constituents of the pair in real 
space, thus rendering them useful as a quantum resource 
of entanglement. While most of the previous analysis of 
entanglement was restricted to the tunneling limitiiSii, 
we have overcome this restriction and have demonstrated 
universal violation of Bis in setups based on electron cor- 



relators. We have determined the degree and duration in 
time of the BI violation and have found pronounced de- 
pendencies on the choice of observables. 
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